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Abstract 



A non-conventional approach to calculating reactions in quantum mechanics 
is presented. Reaction observables are obtained with bound state calculation 
techniques. The accuracy of the method to calculate few-nucleon response 
functions is discussed. 
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I. INTRODUCTION 



The purpose of the method discussed is to solve reaction problems with bound state 
calculation techniques. The approach can be applied to both perturbation and strong inter- 
action induced inclusive and exclusive processes. Below the case of inclusive perturbation 
induced reactions is considered. 

The general procedure was introduced in 1985 The originally proposed Stieltjes trans- 
form led to stability problems 0. In Ref. || it was found that a Lorentzian kernel allowed 
to cure the instabilities. Accurate results on three- and four-nucleon response functions were 
obtained with this kernel @-[7[. 

In Sec. 2 the procedure is outlined. In Sec. 3 the stability issue is explained and the 
inversion technique is described. In Sec. 4 the choice of the transform kernel is discussed. 
In Sec. 5 the accuracy of the method in the "exactly solvable" d(e, e')np problem is demon- 
strated, and its features in the calculations of the three- and four-nucleon response functions 
are discussed. 

II. OUTLINE OF THE METHOD 

We need to obtain the response functions of the form 

R(e) = Jdf\(^ f \d\%)\ 2 S(E f -E -e) (1) 

starting from a few-body Hamiltonian. Here e is the excitation energy, is the ground state 
wave function, tyf is a complete set of the final state wave functions, and O is a transition 
operator. In general the continuum wave functions \I/ f include infinite numbers of various 
channels with three or more fragments. The number of these functions is infinite as well. 
Therefore the calculation of R via direct use of Eq. (1) would in general be impractical. 
Complications may occur even at very low energy due to the Coulomb interaction in channels 
with more than two fragments. The present method circumvents the above problems due 
to the fact that the continuum wave functions do not enter the calculation at all. 

The method proceeds in two steps. At the first step an integral transform of R with a 
smooth kernel K 

$(a) = J K(a,e)R(e)de, (2) 

is calculated instead of R itself. Using sum-rule techniques one can see that 

= (V \&K(cr, H - E )d\V ), (3) 

H and £0 being the Hamiltonian and the ground state energy, respectively. Expression (3) 
is evaluated with bound-state techniques as described below. As a second step Eq. (2) is 
considered as an integral equation to invert and thus one gets R(e) from $(er). 
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III. INVERSION OF THE TRANSFORM AND ISSUE OF STABILITY 



Let us suppose that the transform (3) is evaluated and consider the subsequent inversion 
of Eq. (2). As it is known solutions to such type of equations are unstable with respect to 
high frequency oscillations. Indeed, let us add an oscillating increment AR to the right-hand 
side of Eq. (2). The transform will get the increment A$, and if the oscillations are rapid 
enough A$ may be small even if the oscillation amplitudes are not small. The quantity 
A$ might become less than possible errors in the numerical calculation of the transform 
and the corresponding rapidly oscillating increment AR could not then be discriminated. 
It is known that a regularization procedure can cure the instability. The mathematical 
justification of such procedures can be found e.g. in Ref. ||. The regularization suppresses 
rapid oscillations. The higher an accuracy in $ is the higher might be the frequency level 
at which the suppression occurs and the finer are the details of a true response that can be 
reproduced with an approximate solution. A stability of the solution with respect to changes 
in the regularization parameter would then indicate that the high frequencies suppressed are 
insignificant for reproducing a response. 

The inversion procedure in the calculations [|l]-0] was the following. The response was 
sought for in the form 

N 

R ( e ) = Yl °nXn(e,a) , (4) 

n=l 

where \ n are known functions including nonlinear parameters a. If one substitutes this form 
into the right hand side of Eq. (2) one has 

JV 

= J2 C nXn((T,a) , (5) 
n=l 

where the x n (a,a) are the transformed basis functions. The parameters c n and a are 
obtained fitting the transform of Eq. (5) to the known results obtained by Eq. (3) at 
many a points. The number of functions N plays the role of the regularization parameter 
and is chosen within the above mentioned stability region. 

Features of a response known beforehand such as its low energy behavior and possible 
information on narrow levels could easily be incorporated into the set {Xn}- This increases 
the accuracy of the inversion. The longitudinal (e, e') few-body responses behave in a rather 
peculiar way at small e. This is caused by the fact that lower multipole contributions to 
the continuum part of the responses can have maxima in the threshold region, while higher 
multipoles exhibit maxima only in the quasielastic peak region. We applied the method 
separately to those lower multipole responses and to the sum of all others retained ||J|. 
This substantially increased the accuracy of the results at low energy. The reason is that 
for simple functions \ n one can better describe the behavior of those pieces of the response 
function with sufficiently low N values in Eq. (TJ) than the total response function. 

IV. THE CHOICE OF THE TRANSFORM KERNEL 

At a given accuracy in the transform $ the resulting response R will be most accurate 
if the kernel is chosen as "narrow" as possible. Indeed, only changes in R occurring at the 
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intervals Ae smaller than the range of the kernel may be hard to resolve since they may 
tend to cancel at integrating in Eq. (2). 

The following choice proved to be efficient ||. Let the transform variable a be complex 
and consider the kernel of the Lorentz form 

1 1 

(a* + e) (a + e) [(a R - e) 2 + <jj\ 

a i being "sufficiently small". According to Eq. (3) this leads to the transform 
Eq. (0) can be written as 

= <*(a)|*(a)>, (8) 

where \P is the unique solution to the Schrodinger-like equation with a source, 

(H-Eo-aR + taj)^ = Q, Q = 6^ . (9) 

Since $ and thus (^|^) does exists, ^ is localized! Hence bound-state type methods can 
be applied to solve Eq. @. 

V. CALCULATIONS OF FEW NUCLEON RESPONSES 

The accuracy of the method was tested with the two-body longitudinal d(e, e') response. 
To approach the situation with A=3 and 4, several percent errors were artificially introduced 
into the transform, see ||. Fig. 1 demonstrates that even in presence of the errors the 
response obtained practically coincides with the exact one calculated in the conventional 
way. 

Longitudinal (e, e') response functions ||,|5] and total photoabsorption cross sections |||7J 
were studied for three- and four-nucleon cases. In Q Reid and Bonn realistic NN forces were 
employed. Eq. (|9[) was cast into the form of inhomogeneous Faddeev-type equations. In 
0-0] local NN (Malffiet-Tjon type) potentials reproducing s-wave NN phase shifts up to 
the pion threshold were used. Eq. (|9|) was solved with the help of the Jastrow correlated 
hyperspherical harmonic basis (CHH). The 07 values were taken from the interval 5 MeV< 
o~i < 20 MeV. The smaller 07 the more accurate is the inversion but the larger is the number 
of CHH basis functions required to achieve the same accuracy in the transform. In fact for 07 
tending to zero the scattering regime is recovered. In other problems, such as the response 
of the 11 Li nucleus in the 3-cluster model much lower 07 values can easily be used. For an 
efficient inversion the transform should be obtained in the <7r region covering the interval of 
e values of interest. Besides 07 = const, the choice 07 = o"/(ct_r) is possible. In Refs. 0,0] the 
a 1 values were taken to increase with o~r. This improved the accuracy in R at low energy. 

To judge the quality of the responses obtained the following criteria were applied. First, 
the input transform should be obtained from the dynamic equation ((^) with a sufficient 
accuracy. A typical trend of convergence is demonstrated in Fig. 2. The quantity K max 



4 



governs the number of the CHH basis functions retained at solving Eq. ||). The K max values 
in Fig. 2 are sufficient to provide an accurate response and they are substantially lower than 
those required in the bound-state calculation with the same NN potential. Second, as said 
above, the responses obtained should be stable with respect to a change in iV in Eq. (U) in 
some AiV interval. This is the most important criterion. In the harder 4-nucleon calculations 
we had AiV ~ 3 — 4. In some easier cases the AiV values were much larger. Third, various 
choices of 07 in Eq. (g) and sets {Xn} in Eq. ([|) should lead to approximately the same 
responses provided that the above mentioned criterion is fulfilled, and this also tested in our 
calculations as well. An additional check is provided by sum rules. The quantities / R(e)de 
and / R(e)ede were compared with the sum rule expressions evaluated independently. The 
differences were about 1% or less for the non-energy weighted and less than 2% for the 
energy-weighted sum rules. 

To sum up, the method proved to be a practical tool for studying few-nucleon responses. 
It would be of interest to try the approach also for exclusive and strong interaction induced 
reactions. 
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FIGURES 



FIG. 1. Longitudinal deutron form factor at q 2 = 5 fm 2 . Conventional calculation: full curve; 
from inversion of the Lorentz transform: dashed curve. 

FIG. 2. The Lorentz Transform with various K max values. 
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